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ABSTRACT 

Numerical  aolutlons  of  three -parameter  vortex  flows  over  flat 
surfaces  and  of  rotating  flows  of  von  Karman's  and  Bodewadt's  types 
are  presented  for  a  variety  of  Reynolds  nuiid>ers.  The  results  have 
been  obtained  on  the  basis  of  an  extended  boundary  layer  theory,  which 
allows  a  reduction  of  the  Kavier-Stokes  equations  to  a  set  of  Volterra 
Integral  equations.  The  new  flows  are  In  good  agreement  with  available 
physical  observations. 
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1.  Introduction 

Rotating  flows  normal  to  a  flat  surface  have  been  studied  In  flO] 
and  [11]  on  the  basis  of  an  extension  of  Prandtl's  boundary  layer 
theory.  The  Investigations  In  [10]  dealt  with  the  vortex  motion  of  a 
viscous  fluid,  which  Is  produced  by  a  very  long  rod  normal  to  a  flat 
surface  (see  Figure  1).  In  [11]  the  rotating  flows  of  von  Kirm^n  and 
Bbdewadt,  which  are  generated  by  a  rotating  flat  plate  In  a  fluid  at 
rest  or  by  a  fluid  In  solid-body  rotation  over  a  fixed  plate,  are  re¬ 
examined.  In  all  three  cases  the  Navler-Stokes  equations  have  been 
reduced  to  sets  of  nonlinear  ordinary  differential  equations,  which 
are  connected  with  appropriate  boundary  data. 

In  the  present  paper  the  nonlinear  boundary  value  problems  are 
transformed  Into  equivalent  sets  of  nonlinear  Volterra  Integral 
equations,  which  can  be  solved  by  efficient  Iteration  procedures. 
Complete  numerical  results  for  all  three  problems  have  been  computed 
and  are  displayed  In  the  following  sections  for  a  selected  variety 
of  Reynolds  numbers.  Characteristic  properties  of  these  flows  are 
pointed  out  and  compared  with  actually  observed  phenomena.  In 
particular,  the  properties  of  vortex  flows  over  flat  surfaces  are 
compared  with  phenomena  observed  In  hurricanes.  Although  the  vortex 
models  considered  are  only  very  rough  approximations  to  real  hurricanes, 
the  qualitative  agreement  Is  satisfactory. 
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It  Is  a  distinctive  feature  of  all  flows  considered,  that  they 
depend  on  three  characteristic  parameters.  This  Is  a  result  of  the 
complete  boundary  data  which  have  been  deduced  from  physical  flow 
models  In  order  to  specify  a  unique  solution  of  the  elliptic  Navler- 
Stokes  equations.  Because  of  the  great  Importance  of  this  procedure 
for  the  theory  of  real  flows.  It  seems  worthwhile  to  Insert  here  a 
simple  example,  ;dilch  displays  the  fundamental  dependence  properties 
of  the  solutions  of  the  Navler-Stokes  equations. 

As  Is  shown,  for  Instance,  In  [9]  the  Couette  flows  between  a 
fixed  plate  and  a  moving  parallel  plate  depend  on  the  kinematic  vis¬ 
cosity  of  the  fluid,  the  speed  of  the  moving  plate,  the  distance  of 
the  parallel  plate,  and  the  constant  pressure  gradient.  In  addition 
to  these  four  characteristic  parameters  the  Couette  flows  are  also 
specified  by  the  well-known  similarity  assumption.  Vlhlle  the  simi¬ 
larity  assumption  restricts  the  type  of  the  singularities  which  are 
permitted  at  the  entrance  and  the  exit  of  the  flow  field,  the  con¬ 
stant  pressure  gradient  specifies  their  strengths.  If  the  moving 
plate  Is  removed,  one  obtains  similar  flows  along  an  Infinite  plate 
which  depend  on  three  characteristic  parameters. 

It  Is  one  of  the  most  significant  features  of  the  extended 
boundary  layer  theory  and  the  Integration  procedure  developed  In 
[10]  that  the  correct  dependence  properties  of  real  flows  are  pre¬ 
served.  The  classical  boundary  layer  theory  loses  this  property  by 
truncating  the  elliptic  Navler-Stokes  equations  to  parabolic 
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differential  equatlone.  This  may  be  demonetrated  by  a  simple  example,  | 

i 

which  shows  the  basically  different  dependence  properties  of  solutions 
of  elliptic  and  parabolic  partial  differential  equations. 

The  elliptic  Laplace  equation 

U^  +  Uyy  -  0  (*) 

yields  regular  solutions  In  the  unit  square  [O^x^l,  0^y£l]  for 
any  boundary  data,  which  are,  e.g.,  piecewise  continuous.  However,  the 
parabolic  differential  equation 

yy  ^ 

with  the  general  solution 

u  -  f(x)y  +  g(x) 

does  not  allow  the  prescription  of  boundary  values  along  x  ■■  0  and 
X  -  1  In  addition  to  data  at  y  >  0  and  y  ■  1.  For  instance, If  the 
boundary  data  u(x,0)  -  u(x,l)  *  0,  u(0,y)  ■  (y),  and  u(l,y)  -  GgCy) 

are  Imposed,  no  solution  exists  to  the  parabolic  equation  (**}  no 
matter  how  small  |ei(y)|  4  0  and  lea(y)|^  0  may  be  assumed. 

In  the  problems  of  von  Kirman  and  Bbdewadt  uniqueness  has  been 
achieved  by  very  restrictive  similarity  assumptions,  which  Indirectly 
determine  the  well-known  sink  and  source  distributions  of  Infinite 
Initial  velocities  at  large  llstances  from  the  axes  of  rotation. 

Consequently,  the  physical  models  of  these  flows  remain  unknown. 
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This  has  first  been  pointed  out  by  Stewsrtson  (see  [4,  5,  13])  idio 
questioned  the  physical  meaning  of  the  oscillating  Bbdewadt  flow. 

Moore  (see  [4,  5])  doubted  also  the  usefulness  of  the  Bbdewadt 
so'lutlon,  because  the  flow  shows  no  tendency  of  boundary  layer 
separation  from  the  surface  at  the  axis  of  rotation.  This  should  be 
expected  as  the  secondary  flow  Is  of  wake  type  (see  [12]).  Indeed, 
simple  experiments  In  a  tea  cup  (see  [4])  display  very  clearly  a 
separation  of  the  fluid  from  the  bottom  of  the  cup.  When  the  stirring 
of  the  tea  Is  sufficiently  slow,  tea  leaves  move  toward  the  axis  of 
rotation  and  heap  up  at  the  bottom  of  the  cup.  However,  when  the 
stirring  Is  sufficiently  fast,  the  tea  leaves  settle  distinctly  on  a 
ring  around  the  axis  of  rotation.  This  Interesting  phenomenon  Is  a 
graphic  demonstration  of  a  separated  wake  flow  along  the  axis  of  ro¬ 
tation.  It  Is  comparable  with  the  settling  of  dust  around  comers 
of  rooms  which  are  not  well  kept.  While  the  dust  settles  on  rings 
around  the  corners,  the  comers  themselves  remain  clean.  This  also 
Indicates  a  separated  wake  In  the  comer  Into  which  no  dust  particles 
can  enter. 

Another  peculiar  result,  which  concerns  both  the  von  KiLrm^n  and 
the  Bbdewadt  flows.  Is  the  constant  "deviation  angle"  of  the  spiral 
motion  from  a  circular  motion  at  the  surface  of  the  plate  (see  [9]). 
Since  the  secondary  radial  velocity  depends  strongly  on  the  primary 
tangential  velocity  (compare  [10]),  It  should  be  expected  that  the 
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deviation  angle  of  the  atreafflllnea  at  the  plate  would  depend  on  a 
Reynolds  number.  The  flows  of  von  Xirmiln  and  Bodewadt  do  not  ful¬ 
fill  this  expectation,  which  seems  to  indicate  that  these  flows  are 
considerably  affected  by  the  strong  sinks  and  sources  assitmed  at 
large  distances  from  the  axes  of  rotation. 

2.  Vortex  Flows  Normal  to  a  Plate 

An  axlsysnetric  vortex  flow  over  a  flat  surface  has  been  defined  In 
[10]  as  a  solution  of  the  Havier-Stokas  equations 


uu. 


+  -  r  ■  ■  p  fr  *  t  (“),  +  u„  j 

(?)p  * 


UV,  +  WVj  +  ^ 


UWy  +  ww^ 


+  V 


(ru)y  +  (rw)j  -  0 


by  the  boundary  data 


Vrr  + 


w„+-Wr  +w„ 


2«o}: 

z<r|-  u-0,  v-O,  w-0 

•'  TA  log  jr 


-  1 


-0.^-1. -^-1. 
J  TA  log  T* 


(1) 

(2) 

(3) 

(^) 

(5) 

(6) 

(7) 
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In  these  equations  u,  v,  and  w  denote  the  velocity  components  of  the 
vortex  flow  which  correspond  to  the  coaxial  cylindrical  coordinate 
system  (r,  cp,  z).  The  variable  pressure  and  the  constant  density  and 
kinematic  viscosity  of  the  fluid  are  designated  by  p,  p,  and  u.  While 
the  vortex  strength  T  and  the  radial  extension  r^  of  the  axial  loga¬ 
rithmic  sink  (see  Fig.  1)  are  at  one’s  disposal,  the  constant  value  of 
A  must  be  determined  simultaneously  with  the  solution. 

Guided  by  an  extension  of  Prandtl's  boundary  layer  theory  a 
first  order  reduction  of  the  Navler-Stokes  equations  to  an  ordinary 
boundary  value  problem  has  been  found  In  [10],id)lch  Is  valid  In  the 
vicinity  of  the  line  r  >  tQ.  This  approximation  was  achieved  with 
the  aid  of  the  limiting  line  z  *  6(r)  of  the  boundary  layer  along 
the  surface  at  z  >  0.  This  limiting  line  was  defined  by  the  s  - 
condition 

v(r,z)  -  £  (1  -  5)  (9) 

and  remained  to  be  found  together  with  the  primary  tangential  veloc¬ 
ity  V. 

After  applying  the  following  similarity  transformation 

r  “  C  ■  ^  "  Reynolds  number)  (10) 

«  ■  y  U(C).  V  -  £  V(C),  w  -  FA  log  ^  W(0,  f  “  -  2^  I’(C)  (11) 

U  -  -  G(C),  W  -  i  [g(0  -  CG(0]  (12) 
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one  arrives  at  the  following  ordinary  differential  equations 
G  +  o*C6  -  ^  (G*  +  V*  -  P) 

V  +  o*CV  -  0 

G(G  -  CG)  -  — 

2cr® 

Which  must  be  Integrated  under  the  boundary  conditions 

C  -  0:  G  -  0,  G  -  0,  V  -  0, 

C  -  G  -  G„  ,  V  -  1,  P  -  1, 


(13) 

(14) 

(15) 

(16) 
(17) 


Simultaneously  with  this  reduction  one  obtains  the  equation  of  the 
limiting  line  of  the  boundary  layer  In  the  vicinity  of  r  >  ro  In  the 
form  (see  Flg«  1) 


The  constant  value  A  Is 

A  -  2  ^  G„  , 

*0 

idiere  the  characteristic  parameter  a  Is  determined  by  the  e  >  con¬ 
dition 


(19) 


C  -  Is  V  -  1  -  c. 


(20) 
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For  an  Integration  of  the  remaining  ordinary  boundary  value  problem 
it  is  convenient  to  introduce  the  following  new  variablea 

T)  -  oC  g  -  OG, 

This  substitution  leads  to  the  differential  equations 


•jf  +  Tig  -  I  (^  +  V*  -  P) 
i<g  -  T14)  -  7  P 


(22) 

(23) 


with 


VCT)  .  ^  f 


T\ 

7^ 


(24) 


which  must  be  integrated  under  the  boundary  conditions 
'  T)  -  Oj  g  *  0»  i  *  0, 

T]  -  g  -  g»,  P  -  1. 


(25) 

(26) 


The  ordinary  boundary  value  prpblem  defined  by  the  equations  (22) 
through  (26)  has  been  solved  for  the  Reynolds  number  R  ■  10  by  the 
Runge-Kutta  method  tdiich  started  the  integration  with  an  assumed  set 
of  initial  data  that  had  to  be  improved  successively  until  the  boundary 
condition  (26)  was  sufficiently  met.  The  result  of  this  integration 
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has  been  displayed  In  [10],  The  inefficiency  of  this  method  Is  obvious 
because  two  Initial  conditions  must  be  found  vdilch  correspond  to  the 
two  Ignored  boundary  conditions  (26) . 

For  an  efficient  Integration  of  the  remaining  ordinary  boundary 
value  problem  It  Is  convenient  to  transform  the  differential  equations 
(22)  and  (23)  with  the  boundary  conditions  (25)  and  (26)  Into  an 
equivalent  Integral  equation  (see  [14]).  Assuming  that  the  auxiliary 
function 

D(T1)  -  P(T|)  -  I^T])  -  v®(Tl) 

Is  known,  then  the  differential  equation  (22)  reduces  to  the  linear 
equation 

*r+  TIK  -  -  ^  DCIl) 


(27) 


(28) 


which  Is  Integrable  by  quadrature.  Indeed,  with  the  aid  of  the  error 
function  (24)  one  finds  the  general  solution  of  (28)  In  the  form 

gOl)  -  XV(Tl)  +  X  .  j^vCIl)  -  V(t)J  dt  . 


(29) 


The  boundary  conditions  (25)  and  (26)  lead  then  to  the  nonlinear 
Volterra  Integral  equation 


g(Tl)  -  XV(Tl) 


-fj 


(30) 
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for  gCH),  where 


(31) 


The  functions  g(1))  and  PC'!])  are  obtained  by  direct  Integration  of  the 
equations  (30)  and  (23)  under  the  boundary  conditions  (25)  and  (26). 
They  are 


g(Tl) 


g(t)dt 


0 


(32) 


and 


P(Tl)  -  +  g^(Tl) 


2 


T\ 

t4®(t)dt 
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(33) 


where 

CD 

-  P(0)  -  1  -  +  2 1  t  S®(t)dt. 

This  concludes  the  transformation  of  the  boundary  value  problem  under 
consideration  Into  a  set  of  Volterra  Integral  equations.  It  Is,  vice 
versa,  not  difficult  to  show  that  any  solution  of  the  Integral  equations 
Is  a  solution  of  the  boundary  value  problem.  Thus  the  transformation 
applied  Is  an  Identical  transformation. 
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The  Integral  equationa  (30),  (32),  and  (33)  have  been  numerically 
Integrated  for  a  variety  of  Reynolds  numbers.  Some  selected  solutions 
will  be  discussed  In  the  following  section.  All  solutions  have  been 
constructed  by  an  Iteration  procedure  (see  [10]) ,  which  successively 
Improved  an  appropriate  first  approximation.  It  was  found  that  the 
method  Is  very  efficient,  provided  the  first  approximation  Is  suf¬ 
ficiently  close  to  the  correct  solution. 

For  small  Reynolds  xuimbers  the  Iterations  can  be  started  with 

g(Tl)  =  0,  P(T1)  =  1  (3 

as  adequate  first  approximations.  For  larger  Reynolds  nund>ers  the 
Iterations  may  be  started  by  solutions,  which  are  obtained  for  smaller 
Reynolds  numbers.  In  order  to  Improve  the  rate  of  convergence  of  the 
Iteration  procedure,  especially  In  cases  where  the  Iteration  Is 
started  with  a  very  crude  approximation.  It  Is  helpful  to  average  the 
outcomlng  solution  with  the  entering  approximation  of  the  Iteration 
by  means  of  appropriate  weighting  factors. 

3 .  Properties  of  Vortex  Flows 

Axlsymnetrlc  vortex  flows  normal  to  a  flat  surface  have  been  computed 
and  tabulated  for  various  Reynolds  numbers  (see  Tables  1,  2  and  3). 
Examples  have  been  selected  and  plotted  In  the  figures  3  ,  4,  and  1 
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The  ntnuerlcal  reeulta  confirm  the  i^nomene  of  vortex  flows  which 
were  pointed  out  In  [101.  It  may  be  emphasised  that  the  numerical  cal¬ 
culations  of  all  vortex  flows  computed  Indicated  no  symptoms  of  In¬ 
stability  of  the  flow  near  the  line  r  «  ro  (see  Fig.  1  and  compare 
ri2]).  Accordingly,  vortex  flows  seem  to  remain  laminar  and  attached 
to  the  surface  for  all  Reynolds  numbers  at  least  In  the  vicinity  of 
r  ■  rg.  Nevertheless,  since  the  secondary  flow  is  of  "wake  type" 

(see [12])  within  the  cylinder  r  ■  ro,  separation  of  the  flow  should 
be  expected  to  occur  In  the  neighborhood  of  the  vortex  axis.  However, 
outside  Che  cylinder  r  ro  the  secondary  flow  is  of  "stagnation 
type"  and  tends  to  prevent  any  flow  separation.  Thus,  Che  Instability 
In  the  motion,  which  Is  caused  by  a  flow  separation  at  the  axis,  seems 
Co  fade  away  when  the  flow  changes  Its  character.  This  explains 
the  fact  chat  the  radial  velocity  U  remains  free  of  Inflection  points 
which  characterise  unstable  flows. 

As  was  explained  In  [lOl  vortex  flows  normal  to  a  flat  surface  repre¬ 
sent  approximate  models  of  hurricanes,  provided  exterior  disturbances 
of  Che  vortex  flows  other  chan  those  caused  by  Che  surface  of  the 
earth  are  excluded.  Consequently,  the  numerical  results  may  be  com¬ 
pared  with  phenomena  observed  outside  the  cores  of  hurricanes. 

Despite  the  fact  that  real  hurricanes  are  highly  distorted  by  the 
tremendous  rainfall  Inside  Che  circle  r  -  rg,  by  Che  change  of  the 
density  and  Che  turbulence  of  the  air,  etc.,  the  agreement  with 
available  observations  appears  to  be  satisfactory. 
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lnd««d,  th«  obacrvatlont  confirm  the  wake  and  stagnation  charactar 
of  the  secondary  flow,  idiere  the  rainfall  area  coincides  roughly  with 
the  circular  region  r  ^  r,,.  The  rapidly  Increasing  rainfall  toward 
the  core  of  a  hurricane  (see  [15»p..l30])  Is  Indicated  by  the  log¬ 
arithmic  Increase  of  the  axial  valoclty  toward  the  vortex  axis.  The 
strong  dependence  of  the  secondary  flow  upon  the  Reynolds  nusiber  R 
(  see  Fig.  3.4.  and  5)  Indicates  a  radial  shear  stress  at  the  sur¬ 
face,  which  Is  large  compared  with  the  tangential  shear  stress. 

This  explains  the  very  strong  radial  ocean  waves  produced  by  hurri¬ 
canes  which  are  known  as  ocean  swells  (see  [7 ,  p.298]).  As  cen  be 
deduced  from  figure  3 ,  Inside  the  cylinder  r  •  Tq  the  secondary  flow 
near  the  surface  converges  to  the  surface  before  It  finally  leaves  the 
vicinity  of  the  surface  In  the  normal  direction  (see  Fig.  1).  Accord¬ 
ing  to  Rlehl  this  remarkable  phenomenon  has  also  been  observed  In 
hurricanes  (see  [7  ,  p.  320]) . 

In  order  to  get  a  rough  Idea  of  the  boundary  layer  thickness  near 
the  core  of  a  hurricane,  one  may  consider  a  hurricane  of  average  slse, 
which  has  a  rainfall  araa  of  about  2ro  w  2000  km  In  diameter  (see 
[15, p. 130]).  If  the  limiting  line  of  the  boundary  layer  s  >  d(r) 

(see  Eq.  (18))  Is  determined  to  a  relative  accuracy  of  ■>  1%,  then 
the  characteristic  parameter  o  was  found  to  bo  a  m  2.5  (see  [10]). 

With  these  assumptions  the  boundary  layer  thickness  at  the  edge  of 
the  hurricane  core,whlch  la  at  r  30  km  from  the  vortex  axis  (see 
[  7  p.297]),  can  be  computed  by  equatlon(18),  idilch  yields  6  9  km. 
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This  rsBult  axplains  the  observed  thick  boundary  layers  (see  [6]), 
Which  are  produced  by  the  friction  forces  along  the  surface  of  the 
earth  In  hurricanes. 

4.  Von  Kirtn^n  and  Bbdewadt  Flows 

The  rotating  flows  of  von  K^m^n  (-y  >  0)  and  Bbdewadt  (y  •  1)  have 
been  defined  In  [11]  as  solutions  of  the  Navler-Stokes  equations  (1) 
through  (4)  by  the  boundary  data  (see  Fig.  2). 


|r|  <  00  • 

»  e 

C  ■  0  ^ 


u  ■  0,  V  ■  (1  -  y)*"  r,  w  ■  0 


}:  . 


0,  V  ■  Y  <"  W  “  Woe 


|r|  - 

> : 

z  >0 


U  _  Y 

■.  *  Fr  HTT  Y»  w  -♦  w»  , 


In  both  problems  the  angular  velocity  (u  and  the  constant  Ws>  which 
determines  the  strengths  of  the  singularities  at  the  points  (r  ■  ±  oo, 
z  ■  0),  are  at  one's  disposal. 

After  Introducing  the  limiting  line  of  the  boundary  layer  In  the 


6(r)  “a  -  br®  +  ...  , 


which  Is  determined  by  the  c-condltlon 


z  -  6(r):  SilaSi  -  y  +  (.  ijY 

n\m  '  w 
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the  Navler-Stokee  equations  can  be  reduced  around  r  ■■  0  to  a  set  of 
ordinary  differential  equations.  This  first  order  reduction  requires 
the  following  similarity  transformation 


r  -  r,  C  - 


ufi 


6(r)  »  P 


u 


•  Wtx>® 

*  V  -  HO 

(41) 

),  w  -  Woo  W(0 

(42) 

(43) 

and  the  following  conventions  abou(  the  constants 
■  7  §  *  o*  “ 

The  similarity  transformation  leads  to  the  ordinary  differential 
equations  (see  [11]) 


(44) 


G  +  2tf®(4C  -  RG)G  +  CJ®R(G®  -  V®  +  y  -  4o®CH)  -  0 
V  +  2o*  (4C  -  RG)V  +  2ci*RGV  -  0 
G  +  2o*(2£  -  RG)G  +  C3®RH  -  0, 

which  must  be  Integrated  under  the  boundary  conditions 


(45) 

(46) 

(47) 


C  -  0:  G  -  0,  G  -  0,  V  -  1  -  Y  (48) 

C  -  G  -  G*,  V  -  Y»  H  -  Goo".  (49) 


Vfhen  the  solution  Is  found  for  a  specified  Reynolds  number  R,  the 
e'condltlon 

C  -  1:  V  -  Y  +  (-  1)"^  e  (50) 
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determines  the  characteristic  parameter  o.  Thus,  tte  unknown  coefficients 
a  and  b  of  the  parabolic  limiting  line  of  the  boundary  layer  (39)  and  the 
corresponding  can  be  computed  from  equations  (44).  As  was  shown  in 
[ll],  in  physical  applications  the  ratio  a/b  may  roughly  be  computed 
by 


where  @  represents  the  diameter  of  an  appropriate  finite  flow  model 
(see  Fig.  2). 

The  corresponding  finite  flow  modelb,  which  were  described  in 
[ll]  and  which  are  sketched  in  figure  2,  show  graphically  the  dependence 
of  the  von  Karman  and  Bodewadt  flows  on  3-parameter  Reynolds  numbers. 
Indeed, if  the  slit  between  the  cylinder  and  the  disk  is  of  negligible 
size,  the  finite  flow  model  depends  on  four  essential  parameters; 
the  kinematic  viscosity  v  of  the  fluid,  the  angular  velocity  uu  of  the 
rotating  disk  or  the  rotating  cylinder,  and  the  height  or  and  the  diam¬ 
eter  3  of  the  cylinder  containing  the  fluid,  vnien  carrying  out  the 
limits  a  and  3  one  parameter,  for  instance,  the  ratio  ar/3 

emerges  as  a  characteristic  parameter  for  the  infinite  flow  model. 

In  accordance  with  the  Reynolds  number  (44)  it  is  useful  to  replace 
the  ratio  a/&  by  the  ratio  a/b  as  a  parameter  at  one's  disposal. 

Hence,  a  finite  flow  modal  may  be  considered  similar  to  an  infinite 
model  if  they  both  agree  in  the  ratios  a/b,  that  is,  in  the  products 
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of  the  boundary  layer  thlcknesa  a  and  the  radius  of  curvature  l/2b 
at  the  axea  of  rotation  (see  equation  (39)). 

As  In  the  vortex  problem  (Section  2)  It  Is  helpful  to  transform 
the  remaining  boundary  value  problem  Into  equivalent  Volterra  Integral 
equations.  This  transformation  may  be  carried  out  after  Introducing 


the  following  new  scales  for  all  variables  concerned: 

"H  “  crC,  g  ■  oG,  h  ■  o**H  •  ^^2) 

In  the  new  variables  the  equations  (45),  (46),  and  (47)  assume  the 
following  form 

g-+  STlg  -  -  RA(T1),  A(T))  -  Y+  ^-2gii"V®-4Th  (53) 

V  +  g-nv  -  -  RB('n),  B(T))  -  2(gV  -  gV)  (54) 

h  -  2gg  -  i  (g  +  4T1  g)  ,  (55) 

These  equations  must  be  integrated  under  the  boundary  conditions 

Tl-0:  g-0,  g-0,  V-l-Y  (56) 

"n  -  8  -  8a>  ~  finite,  V  -  y,  h  •  (57) 

With  the  aid  of  the  error  function 

-2T1 

E(T1)  -  erf(2Tl)  -  ^  I  e-'^^dt  ,  (58) 


the  solutions  of  the  boundary  value  problem  obtained  are  exactly  the 
solutions  of  the  following  set  of  Volterra  Integral  equations: 
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-T1 

g(Tl)  -  XiE(71)  -  R  [E(T1)  -  E(t)]  dt 

Jo  E(t) 

.T1 

veil)  -  1  -  Y  +  XaE(Tl)  -  R  ^  [E(T1)  -  E(t)3dt 

Jo 

h(Tl)  -  X3  +  g2  -  I  -  I  tg(t)dt 

Jo 


(59) 


(60) 


(61) 


with 


I 


Ajc) 


(62) 


2y  -  1  + 


“I 


(63) 


r  “ 

■r| 

•'o 


tg(t)dt 


(64) 


The  solutions  of  the  Integral  equations  (59)  through  (61)  may  be 
obtained  by  an  Iteration  procedure  which  successively  Improves  suit¬ 
able  first  approximations.  Numerical  calculations  have  shown  that 
for  small  Reynolds  numbers  R  the  Iterations  can  be  started  with 

g  s  0,  h  s  0,  V  -  1  -  Y  +  (2y  -  l)E('Il),  (65) 

For  larger  Reynolds  numbers  R  the  Iterations  may  be  Initiated  by 
solutions,  which  are  obtained  for  smaller  Reynolds  numbers. 

5.  Properties  of  von  Karman  and  Bodewadt  Flows 

Rotating  flows  of  von  Karman  and  Bc^ewadt  have  been  computed  and 
are  presented  In  the  tables  4  through  11  and  In  the  figures  6  through  12. 
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The  numerical  calculations  of  von  Karman's  flows  Indicated  no 
tendency  of  flow  separation  near  the  axis  of  rotation  for  any  Reynolds 
number.  Since  the  secondary  flow  is  of  stagnation  type  (see  [l2])all 
velocity  profiles  are  free  of  inflection  points  which  characterize 
unstable  flows  (see  [9,  12]).  Accordingly,  the  boundary  layers  of 
von  Kerman's  type  are  considerably  decreasing  in  thickness  as  the 
Reynolds  number  Increases  (see  Fig.  7). 

In  contrast  to  the  von  Karman  flows  the  fluid  motions  of 
Bodewadt's  type  remain  stable  only  below  a  certain  critical  Reynolds 
number  ,  which  has  been  computed  to  be 

Rg  8.  (67) 

For  Reynolds  numbers  R  ^  8  the  Bodewadt  flows  are  not  oscillating 
and  Indicate  no  tendency  of  boundary  layer  separation  at  the  axis  of 
rotation.  For  Reynolds  numbers  R  ^  9  no  proper  nonoscillating  flows 
exist,  which  indicates  the  existence  of  a  separated  flow  around  the 
axis  of  symmetry  (compare  [l2]).  It  may  be  mentioned  that  this  phenom¬ 
enon  has  been  carefully  checked  by  changing  the  input  data  vdiich 
determine  the  accuracy  of  the  numerical  method  applied.  In  addition, 
the  results  have  been  rechecked  by  the  Runge-Kutta  method  which  has 
been  applied  to  solve  the  equivalent  system  of  differential  equations. 
The  high  accuracy  of  both  numerical  methods  and  the  stability  of  the 
solutions  may  be  displayed  by  presenting  the  corresponding  results  of 
both  methods  in  table  form  (see  Tables  12  and  13)  for  the  Reynolds 
number  R  >■  1.  In  this  connection  it  may  be  mentioned  that  the 
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Runge-Kutta  method  used  only  five  digit  Initial  values,  which  were 
obtained  by  the  Integral  equation  method.  It  can  be  seen  that  even 
the  higher  order  derivatives  of  the  velocity  profiles  show  no  symptoms 
of  Inflection  points  or  even  oscillations.  The  same  has  been  found  for 
the  Reynolds  numbers  5  through  8.  However,  for  the  Reynolds  number 
R  -  9  the  remarkable  stability  of  Bodewadt's  flows  Is  clearly  dis¬ 
continued.  As  explained  In  [12]  this  result  was  anticipated,  because 
the  Bodewadt  flows  are  of  wake  tjrpe,  for  which  boundary  layer  separation 
should  be  expected.  It  Is,  Indeed,  a  justification  of  the  doubts  which 
were  raised  by  Stewartson  and  Moore  against  the  solution  of  Bodewadt 
(see  Section  1).  Furthermore,  this  significant  phenomenon  can  easily 
be  confirmed  by  the  tea -cup  experiment  described  In  Section  1. 

Due  to  the  wake  character  of  the  Bodewadt  flows,  the  boundary  layer 
thickness  Is  roughly  constant  with  respect  to  Increasing  Reynolds  numbers 
(see  Fig.  11).  As  was  shown  In  [lO],  this  property  Is  also  shared  by  the 
vortex  motions  over  flat  surfaces  (see  Fig.  4).  Furthermore,  almost 
Invariant  dimensionless  tangential  velocities  V  can  be  observed  In  both 
the  solid-body  rotation  and  the  vortex  motion  along  flat  surfaces. 

Special  attention  may  be  given  to  the  secondary  radial  velocities 
of  the  flows  of  von  Karman  (Fig.  6)  and  Bodewadt  (Fig.  10)  and  the 
vortex  flows  over  flat  surfaces  (Fig.  3).  In  all  three  cases  agreement 
can  be  seen  In  the  property  that  the  dimensionless  radial  velocities 
U  are  very  rapidly  Increasing  with  growing  Reynolds  nuiid>ers.  Thus, 
while  the  tangential  shear  stress  at  the  surface  Increases  relatively 
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The  numerical  calculations  ol  von  Karman'a  flows  Indicated  no 
tendency  of  flow  separation  near  the  axis  of  rotation  for  any  Reynolds 
number.  Since  the  secondary  flow  Is  of  stagnation  type  (see  [l2])all 
velocity  profiles  are  free  of  Inflection  points  which  characterize 
unstable  flows  (see  [9,  12]).  Accordingly,  the  boundary  layers  of 
von  lU(zWn's  type  are  considerably  decreasing  In  thickness  as  the 
Reynolds  number  Increases  (see  Fig.  7). 

In  contrast  to  the  von  Karman  flows  the  fluid  motions  of 
Bodewadt's  type  remain  stable  only  below  a  certain  critical  Reynolds 
number  ,  which  has  been  computed  to  be 

Rj.  8.  (67) 

For  Reynolds  numbers  R  ^  8  the  Bodewadt  flows  are  not  oscillating 
and  Indicate  no  tendency  of  boundary  layer  separation  at  the  axis  of 
rotation.  For  Reynolds  numbers  R  ^  9  no  proper  nonoscillating  flows 
exist,  which  Indicates  the  existence  of  a  separated  flow  around  the 
axis  of  synmetry  (compare  [l2]).  It  may  be  mentioned  that  this  phenom¬ 
enon  has  been  carefully  checked  by  changing  the  Input  data  vdilch 
determine  the  accuracy  of  the  numerical  method  applied.  In  addition, 
the  results  have  been  rechecked  by  the  Runge-Kutta  method  which  has 
been  applied  to  solve  the  equivalent  system  of  differential  equations. 
The  high  accuracy  of  both  numerical  methods  and  the  stability  of  the 
solutions  may  be  displayed  by  presenting  the  corresponding  results  of 
both  methods  In  table  form  (see  Tables  12  and  13)  for  the  Reynolds 
number  R  *  1.  In  this  connection  It  may  be  mentioned  that  the 
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slowly,  the  radial  shear  stress  grows  relatively  fast.  Figure  13  shows 
the  "angle  of  deviation  x"  between  the  spiral  motions  and  an  undisturbed 
circular  motion  at  the  surface,  which  Is  defined  by 


tan  X 


-  lim 

z  -  o  v(s,z) 


U'(0) 

V'(0) 


(68) 


In  all  three  cases  the  absolute  value  of  x  starts  at  zero  as  the  Reynolds 
number  Increases  from  zero,  which  Is  In  contradiction  with  von  Karman's 
and  Bodewadt's  solutions.  However,  the  new  results  are  physically 
plausible. 

Finally,  an  Interesting  pressure  distribution  across  the  flow  field 
can  be  seen  In  both  the  von  Karman  (Fig.  9)  and  the  Bodewadt  flows 
(Fig.  12).  In  a  pure  stagnation  flow  without  rotation  (see  1^12])  the 
pressure  Is  monotonlcally  Increasing  toward  a  pressure  high  at  the  surface. 
The  same  Is  true  for  von  Karman's  flows,  except  that  near  the  surface  the 
pressure  Is  rapidly  decreasing  In  order  to  attain  a  pressure  low  at  the 
surface.  Similarly,  In  pure  wake  flows  (see  [l2])  the  pressure  Is 
monotonlcally  falling  toward  a  low  at  the  surface.  Bodewadt's  flows 
follow  the  same  pattern  up  to  a  point  near  the  surface,  from  which 
point  on  the  pressure  Is  rapidly  rising  toward  a  high  at  the  surface. 

These  phenomena  appear  also  In  von  Karman's  solution  (see  [9])  and  In 
vortex  flows  over  flat  surfaces  (see  Fig.  5). 
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•lowly,  the  radial  shear  atreee  grows  relatively  fast.  Figure  13  shows 
the  "angle  of  deviation  x"  between  the  spiral  motions  and  an  undisturbed 
circular  motion  at  the  surface,  which  Is  defined  by 


tan  X  -  -  lim 

z  -  o  v(i;,z)  V  (0) 


In  all  three  cases  the  absolute  value  of  x  starts  at  zero  as  the  Reynolds 
number  Increases  from  zero,  which  Is  in  contradiction  with  von  Karmen's 
and  Bodewadt's  solutions.  However,  the  new  results  are  physically 
plausible. 

Finally,  an  Interesting  pressure  distribution  across  the  flow  field 
can  be  seen  In  both  the  von  Karman  (Fig.  9)  and  the  Bodewadt  flows 
(Fig.  12).  In  a  pure  stagnation  flow  without  rotation  (see  1^12])  the 
pressure  Is  monotonlcally  Increasing  toward  a  pressure  high  at  the  surface. 
The  same  Is  true  for  von  Karman 's  flows,  except  that  near  the  surface  the 
pressure  Is  rapidly  decreasing  In  order  to  attain  a  pressure  low  at  the 
surface.  Similarly,  In  pure  wake  flows  (see  [l2])  the  pressure  Is 
monotonlcally  falling  toward  a  low  at  the  surface.  Bodewadt's  flows 
follow  the  same  pattern  up  to  a  point  near  the  surface,  from  which 
point  on  the  pressure  Is  rapidly  rising  toward  a  high  at  the  surface. 

These  phenomena  appear  also  In  von  Kerman's  solution  (see  [9])  and  In 
vortex  flowe  over  flat  surfaces  (see  Fig.  5). 
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APPENDIX  A 


TABLE  1 


RADIAL  VELOCITY  -  U  OF  VORTEX  FLOWS 


TABLE  2 


AXIAL  VELOCITY  W  OF  VORTEX  FLOWS 


roA 

a 

.1923 

1.204 

1.523 

1.767 

1.857 

1.951 

A 

1 

10 

20 

50 

100 

500 

0 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.1 

-.0087 

-.0102 

-.0131 

-.0153 

-.0187 

-.0263 

.2 

-.0312 

-.0355 

-.0433 

-.0471 

-.0519 

-.0456 

.3 

-.0618 

-.0681 

-.0787 

-.0794 

-.0787 

-.0454 

.4 

-.0950 

-.1009 

-.1053 

-.1029 

-.0917 

-.0372 

.5 

-.1258 

-.1284 

-.1280 

-.1134 

-.0900 

-.0238 

.6 

-.1499 

-.1463 

-.1366 

-.1094 

-.0748 

-.0035 

.7 

-.1638 

-.1516 

-.1327 

-.0912 

-.0476 

+.0252 

.8 

-.1651 

-.1429 

-.1128 

-.0599 

-.0102 

.0623 

.9 

-.1525 

-.1197 

-.0776 

-.0169 

+.0360 

.1067 

1.0 

-.1256 

-.0827 

-.0287 

+.0358 

.0894 

.1573 

1.1 

-.0851 

-.0332 

+.0209 

.0965 

.1487 

,2126 

1.2 

-.0322 

+.0269 

.0889 

.1629 

.2122 

.2714 

1.3 

+.0311 

.0952 

.1508 

.2332 

.2784 

.3323 

1.4 

.1024 

.1695 

.2290 

.3054 

.3460 

.3940 

1.5 

.1793 

.2472 

.3091 

.3777 

.4133 

.4554 

1.6 

.2592 

.3261 

.3756 

.4486 

.4793 

.5154 

1.7 

.3398 

.4040 

.4535 

.5167 

.5427 

.5731 

1.8 

.4190 

.4792 

.5276 

.5809 

.6026 

.6278 

1.9 

.4949 

.5502 

.5966 

.6406 

.6584 

.6791 

2.0 

.5661 

.6160 

.6496 

.6950 

.7095 

.7264 

2.1 

.6317 

.6759 

.7073 

.7440 

.7556 

.7696 

2.2 

.6910 

.7294 

.7583 

.7875 

.7966 

.8082 

2.3 

.7436 

.7765 

.7957 

.8255 

.8326 

.8424 

2.4 

.7896 

.8173 

.8346 

.8582 

.8639 

.8721 

2.5 

.8291 

.8522 

.8675 

.8861 

.8905 

.8975 

2.6 

.8626 

.8814 

.8948 

.9094 

.9130 

.9190 

2.7 

.8906 

.9057 

.9138 

.9287 

.9317 

.9369 

2.8 

.9136 

.9256 

.9325 

.9444 

.9471 

.9516 

2.9 

.9323 

.9416 

.9476 

.9570 

.9594 

.9635 

3.0 

.9474 

.9543 

.9576 

.9668 

.9692 

.9729 

3.1 

.9593 

.9643 

.9673 

.9744 

.9767 

.9801 

3.2 

.9686 

.9721 

.9748 

.9802 

.9823 

.9853 

3.3 

.9758 

.9780 

.9806 

.9844 

.9863 

.9887 

3.4 

.9813 

.9825 

.9845 

.9875 

.9890 

.9906 

3.5 

.9854 

.9859 

.9882 

.9896 

.9908 

.9913 

I 


XABLE  3 


PRESSURE  DISTRIBUTION  P  OF  VORTEX  FLOWS 


1 

10 

20 

50 

100 

500 

0 

.9992 

.9585 

.9172 

.8528 

.8027 

.7141 

.1 

.9992 

.9585 

.9170 

.8526 

.8023 

.7130 

.2 

.9992 

.9580 

.9158 

.8504 

.7989 

.7075 

.3 

.9992 

.9563 

.9118 

.8442 

.7909 

.7004 

.4 

.9991 

.9528 

.9058 

.8340 

.7794 

.6942 

.5 

.9990 

.9474 

.8959 

.8210 

.7669 

.6897 

.6 

.9988 

.9405 

.8860 

.8076 

.7558 

.6877 

.7 

.9986 

.9328 

.8739 

.7957 

.7480 

.6890 

.8 

.9984 

.9252 

.8631 

.7872 

.7444 

.6942 

.9 

.9982 

.9186 

.8550 

.7830 

.7457 

.7036 

1.0 

.9980 

.9137 

.8508 

.7838 

.7519 

.7170 

1.1 

.9979 

.9111 

.8505 

.7895 

.7626 

.7339 

1.2 

.9979 

.9109 

.8540 

.7995 

.7771 

.7538 

1.3 

.9979 

.9131 

.8598 

.8131 

.7945 

.7756 

1.4 

.9979 

.9174 

.8696 

.8293 

.8140 

.7987 

1.5 

.9980 

.9234 

.8816 

.8471 

.8346 

.8222 

1.6 

.9982 

.9305 

.8927 

.8655 

.8553 

.8453 

1.7 

.9984 

.9384 

.9066 

.8838 

.8755 

.8674 

1.8 

.9985 

.9464 

.9203 

.9013 

.8946 

.8880 

1.9 

.9987 

.9543 

.9333 

.9175 

.9120 

.9068 

2.0 

.9989 

.9616 

.9433 

.9320 

.9277 

.9235 

2.1 

.9991 

.9683 

.9541 

.9448 

.9414 

.9381 

2.2 

.9993 

.9742 

.9634 

.9558 

.9531 

.9506 

2.3 

.9994 

.9793 

.9701 

.9650 

.9630 

.9611 

2.4 

.9995 

.9836 

.9768 

.9726 

.9711 

.9697 

2.5 

.9996 

.9871 

.9823 

.9787 

.9777 

.9767 

2.6 

.9997 

.9900 

.9866 

.9836 

.9829 

.9823 

2.7 

.9998 

.9923 

.9896 

.9875 

.9870 

.9866 

2.8 

.9998 

.9940 

.9924 

.9905 

.9902 

.9899 

2.9 

.9999 

.9954 

.9945 

.9927 

.9926 

.9924 

3.0 

.9999 

.9965 

.9959 

.9944 

.9944 

.9943 

3.1 

.9999 

.9973 

.9973 

.9957 

.9957 

.9957 

3.2 

1.0000 

.9980 

.9983 

.9967 

.9967 

.9967 

3.3 

.9985 

.9990 

.9975 

.9974 

.9974 

3.4 

.9988 

.9995 

.9980 

.9979 

.9980 

3.5 

.9991 

1.0000 

.9985 

.9984 

.9985 

RADIAL 


.1 

0 

.0000 

.05 

.5114-3 

.10 

.8142-3 

.15 

.9527-3 

.20 

.9710-3 

.25 

.9088-3 

.30 

.7998-3 

.35 

.6702-3 

.40 

.5388-3 

.45 

.4175-3 

.50 

.3129-3 

.55 

.2272-3 

.60 

.1603-3 

.65 

.1099-3 

.70 

.0733-3 

.75 

.0477-3 

.80 

.0302-3 

.85 

.0186-3 

.90 

.0112-3 

.95 

.0066-3 

1.00 

.0038-3 

1.05 

,0021-3 

1.10 

.0012-3 

1.15 

.0006-3 

1.20 

.0003-3 

1.25 

.0002-3 

1.30 

.0001-3 

!  I 


TABLE  4 

'  II  OF  VON  FLOWS 


1  10  50  100 


.0000  .0000  .0000  .0000  1 

.5105-2  .4469-1  .1147  .1464 

.8125-2  .6940-1  .1512  .1657 

.9504-2  .7916-1  .1456  .1360  1 

.9683-2  .7860-1  .1218  .0962  ! 

.9060-2  .7166-1  .0935  .0618  ! 

.7971-2  .6143-1  .0676  .0366 

.6678-2  .5015-1  .0467  .0199 

.5366-2  .  3929-1  .0311  .0098 

.4157-2  .2968-1  .0200  .0043  ] 

.3114-2  .2169-1  .0125  .0015  i 

.2261-2  .1537-1  .0076 

.1594-2  .1059-1  .0045 

.1093-2  .0709-1  .0026 

.0729-2  .0462-1  .0015 

.0474-2  .0294-1  .0008  \ 

.0300-2  .1082-1  .0004  I 

.0185-2  .0110-1  .0002 

.0112-2  .0065-1  .0001  | 

.0066-2  ,0037-1  .0001  I 

.0038-2  .0021-1 

.0021-2  .0012-1 

.0012-2  .0006-1 

.0006-2  .0003-1 

.0003-2  .0002-1 

.0002-2  .0001-1 

.0001-2  .0001-1 


I 


TABLE  5 


TANGENTIAL  VELOCITY  V  OF  VON  KABM^  FLOWS 


V 

.1 

1 

10 

50 

100 

0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

.05 

.8875 

.8873 

.8713 

.7776 

.7006 

.10 

.7773 

.7769 

.7470 

.5833 

.4647 

.15 

.6714 

.6708 

.6306 

.4263 

.2990 

.20 

.5716 

.5710 

.5243 

.3050 

.1889 

.25 

.4795 

.4788 

.4291 

.2143 

.1177 

.30 

.3961 

.3955 

.3458 

.1479 

.0719 

.35 

.3222 

.3215 

.2743 

.1003 

.0425 

.40 

.2579 

.2573 

.2141 

.0668 

.0235 

.45 

.2031 

.2025 

.1644 

.0438 

.0112 

.50 

.1573 

.1568 

.1241 

.0281 

.0035 

.55 

.1198 

.1194 

.0921 

.0178 

.60 

.0897 

.0894 

.0672 

.0110 

.65 

.0660 

.0657 

.0481 

.0067 

.70 

.0477 

.0475 

.0339 

.0040 

.75 

.0339 

.0337 

.0234 

.0023 

.80 

.0237 

.0235 

.0159  ‘ 

.0014 

.85 

.0162 

.0161 

.0106 

.0008 

.90 

.0109 

.0108 

.0070 

.0004 

,95 

.0072 

.0072 

.0045 

.0002 

1.00 

.0047 

.0046 

.0028 

.0001 

1.05 

.0030 

.0030 

.0018 

.0001 

1.10 

.0019 

.0018 

,0011 

1.15 

.0011 

.0011 

.0006 

1.20 

.0007 

.0007 

.0004 

1.25 

.0004 

.0004 

.0002 

1.30 

.0002 

.0002 

.0001 

I 


I 


TABLE  6 


AXIAL  VELOCITY  W  OF  VOH  KXRHfb)  FLOWS 


W® 

au) 

-.7643-3 

-.7618-2 

-.5954-1 

-.8366-1 

-.6987' 

X 

.1 

1 

10 

50 

100 

0 

.0000 

.0000 

.0000 

.0000 

.0000 

.05 

.0360 

.0360 

.0407 

.0785 

.1264 

.10 

.1247 

.1249 

.1389 

.2432 

.3592 

.15 

.2419 

.2421 

.2654 

.4234 

.5780 

.20 

.3688 

.3692 

.3990 

.5842 

.7442 

.25 

.4924 

.4929 

.5259 

.7129 

.8562 

.30 

.6045 

.6050 

.6379 

.8088 

.9255 

.35 

.7008 

.7012 

.7316 

.8766 

.9650 

.40 

.7799 

.7802 

.8066 

.922  6 

.9856 

.45 

.8423 

.8426 

.8643 

.9527 

.9953 

.50 

.8898 

.8901 

.9072 

.9718 

.9992 

.55 

.9250 

.9251 

.9381 

.9836 

.60 

.9501 

.9503 

.9597 

.9907 

.65 

.9676 

.9677 

.9744 

.9948 

.70 

.9795 

.9795 

.9841 

.9972 

.75 

.9873 

.9873 

.9903 

.9985 

.80 

.9923 

.9923 

.9943 

.9992 

.85 

.9954 

.9955 

.9967 

.9996 

.90 

.9974 

.9974 

.9981 

.9998 

.95 

.9985 

.9985 

.9990 

.9999 

1.00 

.9992 

.9992 

.9994 

1.0000 

1.05 

.9996 

.9996 

.9997 

1.10 

.9998 

.9998 

.9998 

1.15 

.9999 

.9999 

.9999 

1.20 

.9999 

.9999 

1.0000 

1.25 

1.0000 

1.0000 

TABLE  7 


PRESSURE  DISTRIBUTION  P  OF  VON  KARM^t  FLOWS 


V 

Tl  \ 

.1 

1 

10 

50 

100 

0 

-.2946+5 

-.2954+3 

-3.594 

-.3793 

-.1768 

.05 

+.5683+4 

+.5773+2 

+1.468 

+.9429 

+1.043 

.10 

.2713+5 

.2730+3 

4.375 

1.437 

1.335 

.15 

.3815+5 

.3836+3 

5.740 

1.537 

1.328 

.20 

.4173f5 

.4194+3 

6.080 

1.480 

1.255 

.25 

.4038+5 

.4057+3 

5.797 

1.379 

1.180 

.30 

.3614+5 

.3631+3 

5.186 

1.279 

1.118 

.35 

.3053+5 

.3069+3 

4.450 

1.196 

1.071 

.40 

.2463+5 

.2476+3 

3.716 

1.133 

1.038 

.45 

.1909f5 

.1921+3 

3.057 

1.087 

1.016 

.50 

.1429f5 

.1440+3 

2.505 

1.056 

1.003 

.55 

.1035f5 

.1046+3 

2.067 

1.034 

.60 

.7275+4 

.7377+2 

1.734 

1.021 

.65 

.4969+4 

.5068+2 

1.491 

1.012 

.70 

.3301+4 

.3395+2 

1.320 

1.007 

.75 

.2135+4 

.2232+2 

1.203 

J..004 

.80 

.1344-H4 

.1442+2 

1.125 

1.002 

.85 

.8244+3 

.9226+1 

1.075 

1,001 

.90 

.4924+3 

.5908+1 

1.044 

1.001 

.95 

.2862+3 

.3848+1 

1.025 

1.000 

1.00 

.1616+3 

.2603+1 

1.014 

1.05 

.8830f2 

.1872+1 

1.007 

1.10 

.4646+2 

.1453+1 

1.004 

1.15 

.2318+2 

.1221+1 

1.002 

1.20 

.1065+2 

.1096+1 

1.001 

1.25 

.4218+1 

.1032+1 

1.000 

1.30 

. 1000+1 

.1000+1 

\ 

I 

I 

1 

I 


I 


TABLE 


8 


RADIAL  VELOCITY 


U  OF  BODEWADT  FLOWS 


V 

.1 

1 

5 

8 

0 

.0000 

.0000 

..0000 

.0000 

.05 

.1165-2 

.1165-1 

.0582 

.0907 

.10 

.2054-2 

.2054-1 

.1025 

.1594 

.15 

.2662-2 

.2663-1 

.1329 

.2059 

.20 

.3004-2 

.3004-1 

.1500 

.2316 

.25 

.3111-2 

.3112-1 

.1554 

.2391 

.30 

.3028-2 

.3029-1 

.1514 

.2320 

.35 

.2804-2 

.2806-1 

.1403 

.2141 

.40 

.2490-2 

.2491-1 

.1247 

.1895 

.45 

.2129-2 

.2130-1 

.1068 

.1615 

.50 

.1759-2 

.1760-1 

.0884 

.1329 

.55 

.1407-2 

.1409-1 

.0709 

.1059 

.60 

.1092-2 

.1094-1 

.0551 

.0818 

.65 

.0824-2 

.0825-1 

.0417 

.0613 

.70 

.0604-2 

.0605-1 

.0306 

.0447 

.75 

.0431-2 

.0432-1 

.0219 

.0317 

.80 

.0300-2 

.0300-1 

.0153 

.0218 

.85 

.0203-2 

.0204-1 

.0104 

.0147 

.90 

.0134-2 

.0135-1 

.0069 

.0095 

.95 

.0087-2 

.0087-1 

.0045 

.0061 

1.00 

.0055-2 

.0055-1 

.0028 

.0037 

1.05 

.0034-2 

.0034-1 

.0017 

.0023 

1.10 

.0020-2 

.0020-1 

.0011 

.0013 

1.15 

.0012-2 

.0012-1 

.0006 

.0008 

1.20 

.0007-2 

.0007-1 

.0004 

.0004 

1.25 

.0004-2 

.0004-1 

.0002 

.0002 

1.30 

.0002-2 

.0002-1 

.0001 

.0001 

I 


TABLE  9 


TANGENTIAL  VELOCITY  V  OF  BODEWADT  FLOWS 


n\ 

•  .1 

1 

5 

8 

0 

.0000 

.0000 

.0000 

.0000 

.05 

.1125 

.1126 

.1173 

.1261 

.10 

.2227 

.2230 

.2321 

.2496 

.15 

.3286 

.3291 

.3422 

.3674 

.20 

.4284 

.4290 

.4453 

.4770 

.25 

.5205 

.5212 

.5399 

.5763 

.30 

.6039 

.6046 

.6246 

.6639 

.35 

.6778 

.6785 

.6987 

.7390 

.40 

.7421 

.7428 

.7623 

.8017 

.45 

.7969 

.7975 

.8156 

.8527 

.50 

.8427 

.8433 

.8593 

.8930 

.55 

.8802 

.8807 

.8944 

.9241 

.60 

.9103 

.9107 

.9221 

.9473 

.65 

.9340 

.9343 

.9435 

.9643 

.70 

.9523 

.9525 

.9596 

.9764 

.75 

.9661 

.9663 

.9716 

.9847 

.80 

.9763 

.9765 

.9804 

.9904 

.85 

.9838 

.9839 

.9867 

.9941 

.90 

.9891 

.9891 

.9911 

.9965 

.95 

.9928 

.9928 

.9941 

.9979 

1.00 

.9953 

.9953 

.9962 

.9988 

1.05 

.9970 

.9970 

.9976 

.9994 

1.10 

.9981 

.9981 

.9985 

.9997 

1.15 

.9989 

.9989 

.9991 

.9998 

1.20 

.9993 

.9993 

.9994 

.9999 

1.25 

.9996 

.9996 

.9997 

1.0000 

1.30 

.9998 

.9998 

.9998 

I 


I 


TABLE  10 


AXIAL  VELOCITY  W  OF  BODEWADT  FLOWS 


w„ 

Su) 

.2953-2 

.2955-1 

1.480 

.2251 

_ 

.1 

1 

5 

8 

.00 

.0000 

.0000 

.0000 

.0000 

1 

.05 

.0205 

.0205 

.0204 

.0209 

: 

.10 

.0758 

.0758 

.0755 

.0773 

! 

.15 

.1564 

.1564 

.1558 

.1592 

i 

.20 

.2531 

.2530 

.2520 

.2571 

1 

.25 

.3572 

.3571 

.3558 

.3622 

i 

.30 

.4616 

.4615 

.4599 

.4673 

i 

.35 

.5607 

.5605 

.5587 

.5667 

.40 

.6505 

.6504 

.6485 

.6565 

1 

.45 

.7288 

.7286 

.7268 

.7346 

i 

.50 

.7946 

.7945 

.7927 

.7999 

\ 

i 

.55 

.8481 

.8480 

.8464 

.8528 

.60 

.8903 

.8902 

.8889 

.8944 

.65 

.9226 

.9226 

.9214 

.9260 

.70 

.9467 

.9466 

.9457 

.9495 

.75 

.9641 

.9640 

.9633 

.9663 

.80 

.9763 

.9763 

.9758 

.9780 

.85 

.9848 

.9848 

.9844 

.9861 

.90 

.9904 

.9904 

.9902 

.9914 

.95 

.9941 

.9941 

.9939 

.9948 

1.00 

.9965 

.9965 

.9964 

.9969 

1.05 

.9980 

.9979 

.9979 

.9982 

1.10 

.9989 

.9988 

.9988 

.9990 

1.15 

.9994 

.9994 

.9994 

.9994 

1.20 

.9997 

.9997 

.9997 

.9997 

1.25 

.9999 

.9999 

.9999 

.9998 

1.30 

1.0000 

1.0000 

1.0000 

.9999 

\ 

5 


TABLE  11 


PRESSURE  DISTRIBUTION  P  OF  BODEWADT  FLOWS 


.1 

1 

5 

8 

0 

+.9398+4 

+.9396+2 

3.762 

+1.525 

.05 

+.4034+4 

+.4038+2 

1.631 

+  .6278 

.10 

-.1593+3 

-.1522+1 

-  .0300 

-  .0648 

.15 

-.3222+4 

-.3211+2 

-1.229 

-  .5501 

.20 

-.5248+4 

-.5231+2 

-1.996 

-  .8388 

.25 

-.6376+4 

-.6352+2 

-2.383 

-  .9529 

.30 

-.6773+4 

-.6741+2 

-2.461 

-  .9239 

.35 

-.6619+4 

-.6579+2 

-2.304 

-  .7886 

.40 

-.6087+4 

-.6037+2 

-1.990 

-  .5849 

.45 

-.5331+4 

-.5272+2 

-1.587 

-  .3469 

,50 

-.4478+4 

-.4412+2 

-1.152 

-  .1025 

.55 

-.3626+4 

-.3551+2 

-  .7269 

+  .1290 

.60 

-.2838+4 

-.2757+2 

-  .3435 

+  .3333 

.65 

-.215H4 

-.2066+2 

-  .0127 

+  .5050 

.70 

-.1583+4 

-.1493+2 

+  .2588 

+  .6431 

.75 

-.1131+4 

-.1038+2 

+  .4726 

+  .7498 

.80 

-.7854+3 

-.6902+1 

+  .6350 

+  .8294 

.85 

-.5262+3 

-.4337+1 

+  .7542 

+  .8868 

.90 

-.3482+3 

-.2504+1 

+  .8390 

+  .9269 

.95 

-.2221+3 

-.1234+1 

+  .8974 

+  .9540 

1.00 

-.1362+3 

-.3809 

+  .9367 

+  .9718 

1.05 

-.8073+2 

+.1762 

+  .9622 

+  .9831 

1.10 

-.4569+2 

+.5304 

+  .9785 

+  .9903 

1.15 

-.2395+2 

+.7489 

+  .9885 

+  .9944 

1.20 

-.1093+2 

+.8804 
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FIGURE  1  ;  Scheme  of  the  secondary  flow  produced  by  a  vortex  motion  normal  to  a  flat  surface 


FIGURE  2:  Scheme  of  the  secondary  flow  produced  by  a  rotating  disk  In 
fluid  at  rest. 


FIGURE  3:  The  dimensionless  velocity  components  V  and  W  vs.  the 
dimensionless  variable  T]  for  vortex  flows. 
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FIGURE  4:  The  dimensionless  velocity  component  FIGURE  5:  The  dimensionless 
U  vs.  the  dimensionless  variable  1)  for  vortex  pressure  P  vs.  the  dimensionless 
finwa.  variable  1)  for  vortex  flows. 
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FIGURE  7: 


The  dlmenelonless  velocity  component  V  vs.  the  dimensionless 
variable  T]  for  von  Karman  flows. 
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FIGURE  9:  The  dlmenslraleee  pressure  F  vs.  Che  dimensionless  variable 
T|  for  von  Karman  flows. 


u 

FIGURE  liO  :  The  dimensionless  velocity  component  U  vs.  the  dimensionless 
variable  T]  for  Bodewadt  flows. 
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FIGURE  12: 


The  dimensionless  pressure  P  vs.  the  dimensionless  variable  T| 
for  Bc^ewadc  flows. 
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